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^ ! Abstract 
> . 

, The form factor and the mean-square radius of the pion are calculated analytically 

*^ I from a parametrized form of a qq wave function. The numerical wave function was 

' obtained previously by solving numerically an eigenvalue equation for the pion in a 

. particular model. The analytical formulas are of more general interest than just be 

valid for the pion and can be generalized to the case with unequal quark masses. 
Two different parametrizations are investigated. Because of the highly relativistic 
problem, noticable deviations from a non-relativistic formula are obtained. 



1 Introduction 



Hadrons are composite particles and therefore have a size. A quantitative 
measure of this size is the mean-square radius whose experimental value for 
the pion (tt^) is [1] V < r'^ > = 0.67 ± 0.02 fm. One determines it by first 
measuring the electro-magnetic form factor F{Q'^) for sufficiently small values 
of the (Feynman-four-) momentum transfer = —{Pe~Pe'Y ) and then taking 
the derivative at sufficiently small Q^, i.e. 



2 . . drm 



< >= -6 



Q2=0 



The form factor can also be calculated. One of the most remarkable simplicities 
of the light-cone formalism is that one can write down exact expressions for 
the electro-magnetic form factors. As was first shown by Drell and Yan [2] , it 
is advantageous to choose a special coordinate frame to compute form factors 
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and other current matrix elements at space-hke photon momentum. In the 
DreU frame [3], the photon's momentum is transverse to the momentum of the 
incident hadron and the incident hadron can be directed along the z direction. 
With such a choice the four- momentum transfer is —q^^q^ = Q"^ = ql, and the 
quark current can neither create pairs nor annihilate the vacuum [4]. The 
space-like form factor for a hadron is just a sum of overlap integrals analogous 
to the corresponding non-relativistic formula [2]: 



n f 

with £ ■ — I ~ "^^^ + the struck quark, 

[ k±i — Xiq±, for all other partons. 

The d\p,n] symbolize the convolution over all momentum and helicity space 
arguments of every ^„ and e/ is the charge of the struck quark. This holds for 
any (composite) hadron and any initial or final spins S, but is particularly sim- 
ple for a spin-zero hadron like a pion. The wave functions = "^qq, "^qqg, ■ ■ ■ 
are the Fock-space projections of the eigenstate which for a meson for exam- 
ple arc |^„,e.«m) = T.ii'^qqixi, k±^, Xi)\qq) + ^qqg{xu k±^, \i)\qqg) + ...)• Their 
computation is the aim of the light-cone approach to the bound-state problem 
in gauge theory [4], by solving Hi^cl"^) — ^^|^), with the eigenvalues 
being the invariant mass-squares of the physical mesons. 

Strictly speaking, the above expression for the form factor contains contribu- 
tions from all Fock space sectors. In this work, we restrict ourselves to the 
lowest Fock space projection consisting of a qq pair. There exists a nonzero 
probability of finding the pion in its valence state, which can be calculated 
[5]. It is known empirically, that the form factor at low has essentially 
monopole structure [1]. The mean-square radius is essentially all the informa- 
tion there is for low Q^. In fact, for the nuclcon also, there are various effective 
three quark light-cone descriptions. In these models, without an explicit form 
of the effective potential in the light-cone Hamiltonian for three-quarks, one 
proceeds with an ansatz for the momentum space wave function. Both expo- 
nential [4] and power law [6,4] falloff of this wave function at large k± have 
been used. The low properties of the nucleon like the proton magnetic 
moment iip and its axial coupling have been investigated in these models 
[4]. It is reasonable to assume that the contributions from the higher Fock 
components will only refine this initial approximation [7] . 

We take the starting expression as, 

F{Q^) = jdx(fk±(eiijj{x, k± + {1- + e2tp{x, k± - xqi_)j ijj{x, k±), (3) 



for the purpose of calculating its theoretical mean-square radius. 
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The relation involves only the Lz = Sz = component of the general ud 
wave function, where ip{x, k±) = k±; fl) is the (normalized) probability 

amplitude for finding the quarks with anti-parallel helicities, particularly for 
finding the u-quark with longitudinal momentum fraction x and transversal 
momentum k±, and the d with 1 — x and —k±. Their respective charges are ei 
and 62, respectively, with ei + 62 = 1. 

The M(i-component of the pion (tt"'") is available in numerical form, since it 
has been computed recently in the tJ,-model [8]. But the three-dimensional 
numerical integration of Eq.(3) and it subsequent derivation with respect to 
is cumbersome and may be numerically inaccurate. The aim of the present 
work is therefore to calculate the mean-square radius < > analytically by a 
suitable parametrization of the numerical wave function ijj{x, k±). The general 
procedure outlined in the subsequent sections is applicable also to more general 
cases. 



2 General considerations 



Usually, one is able to write down an integral equation in the three variables 
X and k± for the wavefunction il'{x,k±) [8]. The solution of such an equa- 
tion is numerically nontrivial, among other reasons, because the longitudinal 
momentum fractions arc limited to < x < 1. It is therefore advantageous 
to substitute the integration variable x by another variable — 00 < /c^ < 00 
which has the same range than either of the two transversal momenta k±. For 
unequal quark masses, the transformation is defined as. 



with Ei^2 — y"^i,2 + A;^^ -I- A;^. In the lowest order approximation, the effective 
quark masses mi^2 = Tni^2- 

For equal quark masses m„ = rud = m, the substitution is so simple that it 
can even be inverted, i.e. 

<Q = l(l+ \ \ ^ kl = {m- + kl)^^^ {^) 

2V ^Jw? + kl + kl| x{l-x) 



Formally, the three integration variables kz and k± look like a conventional 
3- vector p = {kz,k±). If one substitutes, in addition, the unknown function 
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ip{x, k±) by an other unknown function (f(p) according to 

N ( p'^\^ 
iP{x,kA_) = ^{K,kA_) ^===[1 + — ] , (6) 

^Jx{l - X) \ / 



one gets an identical integral equation in the three variables p, which looks 
like an integral equation in usual momentum space. In the tJ,-model [8], the 
integral equation is further simplified and at the end looks very simple indeed, 



4m^ + ip{p) 

1 r d^p' f Aw? ^ 

TT^ J W \ {p — p') ^ 



4 

3^ 



/i^ + {p — p'Y ^ 



The equation was solved numerically in [8] for spherical symmetry Lp(p) = ^p(p) 
and for the parameter values m = 1.16, /i = 3.8, and a = 0.6904, where masses 
(and momenta) are expressed in units of T = 350 MeV. The calculated eigen- 
value for the ground state agrees with the pion mass-squared to a high degree 
of accuracy, and is stable with respect to changes in (renormalization). As 
can be seen from Fig. 1, it behaves like a power law ip{p) ~ (1 + {p/PaY)~'^ ■, 
rather than as anticipated in [3] like a Gaussian ^{p) ~ cxp(— (p/p^)^). As 



shown in Appendix. A, a value of n 

Pa, 



2 is more likely than others and a fit of 



1 



Pa = 1-471, 



(7) 



reproduces the numerical wave function quite well as shown in Fig. 1, with 

albeit a comparatively large Bohr momentum Pa- Expressing the latter in a 
length, the value of \/?>/pa = 0.66 fm (similar to the experimental rms) is no 
numerical co- incidence, but the result of a constraint on the procedure in [8]. 

— * — * 

One should emphasize that only 4>{x, k±) but not (p{p) — ip{x, k±) is a prob- 
ability amplitude, and that the two can differ appreciably from each other 
according to Eq.(6), particularly for a large Bohr momentum ~ m as in 
Eq.(7). Disregarding this proviso, and mistaking ip{p) as a probability ampli- 
tude for finding the particle with relative momentum p, the calculation of the 
rms-radius would be trivial: One takes the Fourier transform of Eq.(7), i.e. 
FT[(p{p)] ~ exp (— Par), and calculates its second moment to be 

2 J dr r"^ cxp {—2par) 3 

^'^ J dr r"^ cxp {—2par) p^ 



A 



10.0 




10.01 



10.0 



Fig. 1. The pion wave function ^{p) 
is plotted versus p/ {1.552k) in an ar- 
bitrary normalization. The filled circles 
indicate the numerical results, the open 
circles the fit function. The diamonds 
denote the pure Coulomb solution. 



Fig. 2. The function x{p) is plotted 
versus p and fitted to the three expo- 
nents K = 2, |,1 (triangle, box,circle). 
To guide the eye, ip is included by the 
dotted line. Note the increased scale as 
compared to the left. 



However, as the variable conjugate to the relative momentum p, the quantity 
r is the relative distance of the particles. It is not the radius- distance from the 
common center-of-mass, with respect to which the rms is usually calculated. 
The latter is r/2 for equal mass particles, thus 



3 



It will be referred to as the 'non-relativistic estimate', since the above con- 
struction holds approximately if k1 <^ and <^ m^, thus a; ~ ^ according 
to Eq.(5). The question we pursue in the present work is thus: How large is 
the discrepancy between the non-relativistic estimate of Eq.(8) and the quasi- 
exact form factor according to Eq.(3) and its behaviour for low according 
to Eq.(l), particularly for large Bohr momenta. 

Before proceeding with the computation of the form factor, a number of no- 
tational definitions will be introduced, in terms of which the final results turn 
out to be simple. Once one has (^{'p) in a parametrized form like Eq.(7) one 
can transform back to the variables x and /cj_. Since = + one can use 
Eq.(5) to get 

i + ^ = i + ^'^^r^ + i^^^-^(^'^x)- (9) 

p^ x[l — X) Am'' x[l — X) 

The dimensionless variable s is introduced conveniently as 

s = ^, (10) 

Pa 
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as well as the isolation of the pure x-dependence by 



X{x) = l + s 



2V^__2j_ 

x{\ — x) 



The combination 1 +p^/m^ is then trivially obtained by s = 1. One can thus 
compute the form factor according to Eq.(3), i.e. 

F{(1) = I dxdkl\ilj{x,k^)\'f{x,kr,qL): (12) 
<r2> = y dxdk]_\ip{x,kA_)\^g{x,k±). (13) 

The function f{x, k±; qj_) and g{x, k±) contain all the difficulty in the problem, 



27r 



f{x,k±;q±) 



eiip{x, + (1 - x)qi) + e^ipix, k_L - xqi) 



g{x,k±) = -6 f{x,k±;q±) 



4L=0 



'2^jj(x,k±) 

-- -6 -^f{x,k±;q±) 



(14) 
(15) 



4L=0 



that is the integration over the angle between k± and gl {k±qi_ — k±qi_ cos0). 
Finally the size parameter is introduced, 



S 



cal 



exp 



0.67 fm 



(16) 



as a dimensionless measure of the size. 



3 The parametrization of the wave function 



Since the quarks move highly relativistically (p^ ^ m^) and one cannot dis- 
regard the factor ^1 -\-p^/m? in Eq.(6). One way to account for that is to 
parametrize directly 

x(p)^^(p)i/^=(i+g) ^ (17) 



with the two adjustable parameters k and Ps- We have performed two fits, for 
two fixed values of and have calculated analytically the size parameter S 
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according to Eq.(16). The results are: 



«; = 2, 5 = 0.5588, = 2.075, 5 = 0.4577, 
= I, s = 0.6796, Ps = 1.707, S = 0.4321. 



One needs to consider essentially the integral 



j dxd'^k±\il){x, k±) 



Inserting ip according to Eqs.(6) and (17) gives with s = m/ps 



1 



7V2 J x(i -x)J "^"^J 2 











2k ■ 



Z{x, kl) 

The integration over is trivial, and the integration over kl_ elementary, since 



/ 



dk} 



4p^a;(l — x) 



Z{x, kl) {n - 1) [X{x)] 



n-l ' 



(18) 



with X[x) defined in Eq.(ll). Introducing 

} 1 



by definition, one remains with 



1 

7V2 



2k -I 



A{s). 



(19) 



Here the matter rests, since A cannot be integrated in closed form. 

For to compute the means-square radius, we first evaluate the function / as 
defined in Eq.(14). In the general case one gets 
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f{x,k±;q±) = ^Z^ j c 



- 1)2 (1 - (1 - X)a) ' 



+ 



x(l — x) 4p2x(l — x) 



7 



2-K 



1+^2 l^Jl!^ 



xq± 



21 



x{l — x) 4:pjx{l — X 



Using the abbreviative coefficient functions depending on x and kl, i.e. 



1 — X ql 



6i = 



X Apl ' 
1 l-x 



4 X ' 



ci = - 



CiCi = 



1 fci 

a; 2ps ' 
1 l-x 



2p: 



X 



(20) 
(21) 



and correspondingly 62, C2, 62 and C2 by exchanging x 1 — x, gives 



ei 



62 



(61 + Ci COS^)*^ (62 + C2 COS0)'* 



The integration over can not be carried out in closed form for general values 
of K. 



3. 1 The integrals for k — 2 



The normalization integral is now 



1 



Air 



pIA{s). 



The contribution to / from the up-quark becomes to leading order in c 



/ b bb — cc b cc 

7 = 6" '^^-^ 6^6^' 



Inserting this into the definition of g in Eq.(15) gives 

37r 



g{x, k±) = 



x{l — x)p 



2 2,X 



x\ 



Here we have taken both the quark and antiquark contributions (see Appendix 
B). 
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Fig. 3. For k = 2, the functions A{s) and 
< > (in units of ^) are plotted ver- 
sus s by the sohd and the dashed-dotted 
Une, respectively. 



Fig. 4. For k 



I , the functions 



A{s) and < > (in units of ^) are 
plotted versus s by the solid and the 
dashed-dotted line, respectively. 



With the elementary /^^-integration of Eq.(18), one ends up with 



< r >= 



dx 



x2 _ 9 R{s) 
x{l - x) [X{x)t ~ 10p2 A{s) • 



The two elementary integrals A and R are complicated but straightforward 
and yield 



A{s) = j dx 
R{s 



^ 1 -8- 10s2 + 3s^ + 36(s)s2(8-4s2 + s4)) 





1 



[X{x)Y 8(s2 - 1)3 

dx x^ -136 + 72s2 - 56s^ + 15s^ 



y x{l-x) [X{x)f 48(s2- 1)4 

, , , 32 - 16s2 + 36s^ - 22/ + 5/ 
+36(s) 



48(s2 - 1)4 

The auxiliary function b{s) is an abbreviation for 



bis) = "^^^""(^^^ (22) 
VS^ — 1 



The functions are plotted in Fig. 3. As seen there, the asymptotic behaviour 
for s — > oo, 

Ais) ^ — , R{s) — , 
lbs 32s 
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0.0 0.5 1.0 1.5 2.0 

Fig. 5. The root-mean-square radius of 
the pion, {i.e. V< >/0.67 fm), is plot- 
ted versus s by the sohd hne for k = 2 
and by the dashed hne for k = |. 
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Fig. 6. In 1/99 is plotted versus Inp (cir- 
cle), and compared with 41np (dashed) 
and 31np (dotted line). 



is practically reached for values as small as s = 2. Note that in the asymptotic 
limit s — > cxD, the mean-square radius reaches the correct non-relativistic value 
< >nr= 3/(4|}g). As a final result, the size 



he 



9R{s) 
0.67fmTm ^ 10A{s 



0.8401 



9R{s) 
\ 10A{s 



is plotted in Fig. 5 for the mass value m = 1.16 as function of s. 



3.2 The case k = | 



According to Eq.(19), the normalization integral is for k 

^ = 2nplA{s). 



In leading order in c, the contribution to the rms from the up-quark becomes 

f__3b 15 cc 
7 ~ ~2b^Yl^' 



Inserting the derivatives from Eq.(21) gives (see Appendix C) 



g(x, h 



Stt 



XI 



9 
Z 



15X- 



X 
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With the elementary /c^-integration of Eq.(18), one ends up with 



9 9^ /■ dx 9 Ris) 



8 J x{l-x)[X{x)f 16p2A(s)' 



according to Eqs.(19) and (13). The elementary integrals are now 



. } , 1 2 + 52 + 6(5)^2 (-4 + ^2) 
Ais) = / dx ^ = tV — -■ 



dx x'^ -20 - + 



x) [X{x)f 8(s2 - 1)3 

-16 + 8s2 + 10s^ + 3s6 



+ h{s) 



8(s2 - If 



The auxiliary function h{s) is the same as in Eq.(22). The functions are plotted 
in Fig. 4. Finally, the size function 



he 



9R{s) 0.8401 



O.GTfmTm ^J 16A(s) m \l&A{s) 



9R{s) 



is plotted in Fig. 5 for the mass value m — 1.16 as function of s. 



4 Summary and conclusion 



We have analytically calculated the mean square radius for the pion in the 
light-cone approach in the valence sector. The contributions from the higher 
Fock space sectors are expected to only refine this initial approximation. We 
have used a parametrized form of the wave function in momentum space which 
had been obtained previously by solving the eigenvalue equation for the ef- 
fective light-cone Hamiltonian. The mean square radius is calculated by first 
calculating the form factor at sufiiciently low and then taking its derivative 
with respect to at = 0. We have investigated two different values of the 
parametrization. k — 2 is found to be a better fit than k = |. We have shown 
that the mean square radius deviates noticably from the non-relativistic esti- 
mate because of the relativistic effects. In the asymptotic region (large m^), 
it approaches the correct non-relativistic value. 
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A The numerical wave function of the pion 



The function ip{p) was cumputed in previous work [9] and is tabulated in 
Table A.l. It behaves much like an inverse power (p{p) ~ (1 + {p/PaY)~'^- The 
value of K is closer to 2 than to |, as demonstrated in Figure 6. Unfortunately, 
the maximum momentum of the Gaussian quadratures does not allow for a 
preciser statement. 



B The function g{x, k±) for k — 2 

Let us denote the contribution from the up-quark to the form factor function 
simply f{x, k±;q±). For k, = 2 one gets: 

^ e b 
f{x, k±; Qi) = d(j) = Z^ne -. 

I (6 + CCOS0)' (62_c2)2 

Table A.l 

The calculated pion wave function 



p 




P 




(t){p) 


0.015739 


10.00000 


0.083241 


9 


928756 


0.205995 


9.563656 


0.386376 


8 


592389 


0.628023 


6.908239 


0.936173 


4 


877784 


1.318156 


3.077257 


1.784203 


1 


801300 


2.348743 


1.011136 


3.032468 





554441 


3.867276 


0.297233 


4.902437 





154397 


6.223861 


0.076201 


7.998102 





034485 


10.62105 


0.013264 


15.52088 





003404 
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Its derivative can be calculated quite in general 



f b bb-cc b{-2b^ + c^) + 3bcc 



f b b^-c^ 



b{b'^ - C2) 



In the limit g± — > (thus c — > 0) it becomes to leading order 



Inserting the derivatives from Eq.(21) and taking the contributions from quark 
and anti- quark gives one gets thus with 



g{x, k±) 



37r 



2 



3X 
2^ 



1 — X 



e-L- 



X 



+ 62 



X 



1 — X 



The expression in the round bracket has a contribution (1 — 2x) which is odd 
under the exchange x 1 — x, 



1 — X X 
ei h 62 



[(ei + e2)x^ + - 2a;)] (d + ea)^^ 



X 



1 — X 



x{l — x) 



x{l — x) 



It will vanish in the final integration over x and can be omitted. This leaves 
one with the final expression 



g{x, k±) 



x{l — x) p' 



(ei + 62) 



2 3X 



X 



(B.l) 



C The function g{x, k±) for k = | 



The form factor function for the up-quark 



TV 



2eZf 



(6 + ccos0)2 p 



E[g], with < 



p= (6-c)^(6 + c), 
2c 
+ c 



can be integrated in closed form also for k = 3/2 and be expressed in terms of 



the complete elliptic integrals of the second kind, E[q\ = d9y 1 — q sm^9. 
Since dE[q\/dq — {E[q\ — K[q\) /2q, the derivative is 

/ q E[q] - K[q] p 

f q mq] P' 
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Both derivatives q and p have a potentially dangerous singularity be, 



q —cb + be 
q {b + e)e ' 



and 



p (36 + c)6- (3c + 6)c 



P 



2(62 - c2) 



To cure the problem, one must expand 



E[q]-K[q] _ q 5q' 3 

E[q] 2 16 ^~ b + e 



c 5 c 

' 1 + 



46 + c 



to a sufficiently high order, which upon insertion yields 



/ 



+c f^^ 5c 



+ c 



2(6 + c)2 
-6 



36 + c 



4(6 + c) ; 2(62 _ c2 



1 + 



5c 



+ 



6 + 3c 



_2(6+c)2 \^ 4(6 + c)y 2(62 -c2) 
The leading terms in the coefficient of c now tend to cancel, i. e. 



f_ 
f 



+c (^^ 5c 



+ 



2(6 + c)2 
cc 



6 + c 



36 + c 
4(6 + c) J ~ 2(62 - c2) 

6 



56 3 

+ -TTT, T + 



8(6 + c)2 2(6 -c) (62 -c2) 



3 6 15 cc 
~2 6^"8'62' 



In the last step only terms were kept which survive in the limit c — > 0. Inserting 
the derivatives from Eq.(21) gives for the complete amplitude 



g{x, k±) = 



8x{l — x)p 



(ei + 62) 



9 15X 



X 



(C.l) 



where the same simplifications have been done as in the previous section. 
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